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Abstract 

The  local  structure  of  several  disordered  carbons,  obtained  from  radial  distribution  function  analysis  of  elastic  neutron  scattering  data, 
consists  of  small  isolated  graphene  fragments.  The  effect  of  residual  hydrogen  which  presumably  saturates  the  dangling  bonds  in  low- 
temperature  pyrolyzed  materials  was  studied  by  inelastic  neutron  vibrational  spectroscopy  and  computer  modelling  of  lithiated  polyaromatic 
hydrocarbons.  We  discuss  the  possibility  of  exploiting  the  high-pressure  crystal  phase  LiC2  in  primary  batteries,  by  demonstrating  that  trace 
boron  substitition  greatly  enhances  the  ( meta)  stability  of  this  phase  at  atmospheric  pressure.  ©  1997  Elsevier  Science  S.A. 
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1.  Introduction 

Li-carbon  compounds  overcome  many  of  the  disadvan¬ 
tages  encountered  with  metallic  Li  anodes.  The  large  and 
growing  variety  of  carbon-based  materials  suggests  the  pos¬ 
sibility  of  increasing  the  maximum  Li/carbon  ratio,  in  par¬ 
ticular  by  using  disordered  carbons.  Here  we  summarize  our 
results  on:  ( i )  neutron  diffraction  and  radial  distribution  func¬ 
tion  analysis  to  determine  the  local  structure  of  disordered 
carbons;  ( ii)  computer  simulations  [  1  ]  and  vibrational  spec¬ 
troscopy  directed  towards  elucidating  the  role  of  residual 
hydrogen  in  low-temperature  pyrolyzed  carbons,  and  (iii) 
prospects  for  stabilizing  super-dense  high  pressure-synthe¬ 
sized  LiC2  for  use  in  primary  batteries  [2]. 

2.  Local  structure  of  disordered  carbons 

Extensive  work  by  Dahn  et  al.  [3]  shows  that  non-crys¬ 
talline  carbons  exhibit  even  greater  Li  capacities  than  crystal 
graphite.  As  a  first  step  towards  a  microscopic  understanding 
of  Li  uptake  in  such  materials,  we  measured  the  local  struc¬ 
tures  of  a  variety  of  disordered  carbons  using  neutron  radial 
distribution  function  (RDF)  analysis  [4,5]  and  real  space 
modeling.  Representative  data  from  two  of  Dahn’ s  materials 

*  Corresponding  author. 

0378-7753/97/$  1 7.00  ©  1997  Elsevier  Science  S.A  All  rights  reserved 
P//S0378-775  3t  96  )  025  63-3 


are  presented;  their  RDF’s  turn  out  to  be  typical  of  mea¬ 
surements  from  a  wide  variety  of  disordered  carbons, 
obtained  from  different  precursors  with  different  pyrolysis 
temperatures. 

In  RDF  analysis,  the  differential  correlation  function  D{  r) 
is  obtained  by  Fourier-transforming  the  elastic  structure  fac¬ 
tor  S{Q) 

•jC 

D{r)=—  \q[S(Q)-\]  sin(gr)  dQ  =  4irr(p(r)  -p0) 

7TJ 

0 

where  p(  r )  is  the  number  of  atoms  per  unit  volume  at  r,  p0 
is  the  average  microscopic  density  and  Q  =  4tt  sin  0/A  is 
the  momentum  transfer.  Other  quantities  of  interest  are 
the  total  correlation  function  T(r)  =47rrp{  r)  and  the 
RDF  =  4  m^plr),  the  latter  being  the  number  density  at  a 
distance  r  from  a  reference  atom.  For  a  gas,  p(r)  =p0  at  all 
r,  thus  D(r)=  0  and  the  RDF  is  smooth  and  featureless. 
Short-range  (or  long-range)  correlations  in  condensed  matter 
produce  peaks  in  D(  r) ,  and  oscillations  in  the  RDF  about  the 
average  parabola  defined  by  the  r  factor.  Spatial  resolution 
is  limited  by  truncation  broadening  due  to  termination  of  the 
integral  at  the  maximum  experimentally  accessible  momen¬ 
tum  transfer  Qmax.  In  X-ray  diffraction  ( XRD)  0max  is  limited 

o  _  1 

to  ~  1 8  A  using  Mo  Kw,  while  electron  diffraction  is  com¬ 
plicated  by  absorption  and  multiple-scattering.  We  used 
pulsed  dme-of-flight  neutron  diffraction  and  the  glass-liquid- 
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Fig  1.  Neutron  RDF  ( 47 rrp)  of  ( - )  graphite  and  ( - )  LiC^,  obtained 

from  pulsed  time-of- flight  neutron  diffraction 


amorphous  material  diffractometer  ( GLAD )  at  the  intense 
pulse  neutron  source  (IPNS)  of  Argonne  National  Labora¬ 
tory  with  a  (2mu*  =  45  A  ~~  1  [  5  ] . 

Crystalline  graphite  and  LiC6  powders  were  studied  as 
standards.  In  Fig.  1  we  plot  the  RDF's  of  graphite  (dashed 
curve)  and  LiC(1  (solid  curve).  Sharp  peaks  can  be  seen  in 
both,  corresponding  to  shells  of  atoms  at  well-defined  dis¬ 
tances  from  a  reference  atom.  For  graphite,  the  first  peak  is 
at  1.43  A,  corresponding  to  the  smallest  C-C  distance  (in¬ 
plane  bond  length).  The  first  and  second  peak  positions  give 
the  bond  angle,  120±  1°.  The  third,  the  intra-hexagon  dis- 

D 

tance,  is  2.86  A,  twice  the  bond  length,  a  signature  of  planar 
hexagons.  Integrated  intensities  for  the  first  three  peak  are 
2.98.  6.09  and  3.5  atoms,  respectively,  in  reasonable  agree¬ 
ment  with  the  expected  values  3.0,  6.0  and  3.0.  The  first 

o 

interlayer  C-C  correlation  should  occur  at  3.35  A  with  a 
coordination  of  1,  since  only  half  the  C-C  pairs  are  eclipsed. 
This  can  be  identified  as  a  shoulder  on  the  lower  side  of  the 

o 

fourth  ( in-plane )  neighbor  peak  at  3.75  A.  The  average  den¬ 
sity  p().  fixed  by  requiring  the  RDF  to  vanish  below  the  first 
peak,  is  2.26  g/cm\  in  excellent  agreement  with  the  hand¬ 
book  value  2.23  g/cm\  All  of  these  local  observations  are 
entirely  consistent  with  the  well-known  edge-connected 
hexagonal  network  of  sp2-hybridized  carbons. 

Most  of  the  RDF  peaks  in  LiC6  (solid  curve  Fig.  1 )  are 
up-shifted  relative  to  graphite.  This  is  due  to  the  elongation 
of  the  near-neighbor  C-C  distance  upon  intercalation,  a  con¬ 
sequence  of  electron  donation  from  Li  to  carbon  layers  — 
charge  added  to  the  carbon  it  system  slightly  weakens  the  a 

bonds.  Spatial  correlations  between  heteroatoms  Li  and  C 

° 

give  rise  to  the  negative  peak  at  2.36  A  because  the  neutron 
scattering  length  of  Li  (natural  abundance  of  6Li  and  7Li)  is 
negative.  The  observed  Li-C  distance  is  consistent  with  inter¬ 
calated  interstitial  Li  sandwiched  symmetrically  between  two 


eclipsed  hexagons,  which  in  turn  is  consistent  with  the  known 
change  in  graphite  layer  stacking  from  ABAB  to  AaAa, 
where  A  and  B  represent  carbon  layers  and  a:  the  lithium  layer 
[6] .  Similarly,  the  interlayer  C-C  shoulder  at  3.35  A  is  gone 
in  LiC6;  all  the  carbons  are  now  eclipsed  and  the  interlayer 
spacing  has  increased  to  3.7  A. 

The  above  results  on  crystalline  prototypes  give  us  confi¬ 
dence  that  the  corresponding  local  information  will  be  obtain¬ 
able  from  disordered  carbons  using  the  same  methods.  In 
Fig.  2  we  plot  D{  r)  which  emphasizes  the  oscillations  arising 
from  short-range  correlations,  for  two  epoxy  novlac  resins 
( ENR)  pyrolyzed  at  650  or  1000  °C.  Details  of  the  prepara¬ 
tion  can  be  found  in  Dahn's  papers,  see  Ref.  |  3].  These  are 
denoted  ENR650  (H/C  =  0.24,  solid  curve)  and  ENR1000 
(H/C  =  0.04,  dotted  curve).  The  RDF  of  polycrystalline 
graphite  is  shown  for  comparison  ( upper  curve).  The  D(  r)'s 
of  ENR650  and  ENR  1000  are  quite  similar,  with  most  peaks 
located  at  the  same  positions  as  in  graphite,  suggesting  the 
same  bond  length  and  angle,  a  hexagonal  planar  local  motif 
and  sp2-hybridized  carbons.  The  peak  intensities  are  weaker 
than  in  graphite,  the  discrepancy  increasing  with  increasing 
r.  For  the  low-order  peaks,  this  has  a  very  simple  interpreta¬ 
tion  —  a  significant  number  of  edge  atoms  with  only  two 
neighbors.  The  absence  of  strong  peaks  above  ~  6  A  sets  the 
approximate  scale  of  the  in-plane  particle  size.  In  Fig.  3  we 
compare  D(r)  of  ENR  1000  ( solid  curve )  with  a  model  cal¬ 
culation  assuming  isolated  graphene  fragments  5  hexagons 
on  a  side  ( dotted  curve ) .  The  agreement  is  quite  satisfactory. 

We  conclude  that  carbons  obtained  by  pyrolysis  of  epoxy 
novolac  resin  at  temperatures  from  650  to  1000  °C  consist  of 
small  graphene  fragments  with  lateral  dimensions  of  10  A. 
The  fact  that  these  fragments  are  planar  and  the  carbon  atoms 
are  sp2-hybridized  imply  that  the  residual  hydrogens  are 
bonded  to  the  edge  atoms.  Inter-graphene  correlations  are  not 
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Fig  2  Differential  correlation  functions  Di  r  I  of  epoxy  novolac  resm  pyro¬ 
lyzed  at  ( - )  650  °C  ENR650  and  (  ■  •  )  1 000  °C  ENR  1 000.  and  graph¬ 

ite  ( upper  curve  I  for  comparison. 
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Fig  3.  Comparison  between  the  calculated  D(r)  of  {  •  •  •  )  an  isolated 

graphene  fragment  and  ( - )  the  data  from  ENR1000.  The  graphene 

particle  is  shown  in  the  inset 


apparent  in  the  RDF,  indicating  that  the  layer  stacking  ( which 
is  evident  from  non-vanishing  (002)  Bragg  intensity)  is  ran¬ 
dom  parallel  to  the  graphene  fragments,  an  extreme  example 
of  turbostratic  disorder. 


3.  Hydrogen  effect 

Amorphous  carbons  obtained  by  pyrolyzing  organic  solids 
at  T  <  700  °C  retain  substantial  residual  hydrogen  and  exhibit 
surprisingly  large  Li  capacities  [3,7,8].  Furthermore,  an 
empirical  correlation  between  H  content  and  ‘excess’  Li 
capacity  has  been  established  [3].  Using  inelastic  neutron 
scattering  vibrational  spectroscopy,  we  found  that  the  H  sat¬ 
urates  dangling  bonds  on  edge  C  of  the  finite-size  graphene 
flakes  described  above,  similar  to  what  occurs  in  hydrogen¬ 
ated  amorphous  Si.  A  single  H  bonded  to  an  otherwise  two- 
coordinated  edge  carbon  preserves  the  sp2p7T  hybridization, 
and  thus  the  possibility  of  binding  either  a  second  H  to  give 
a  methyl  moiety,  or  Li  which,  if  occurring  on  neighbor  pairs 
of  edge  C,  would  be  analogous  to  the  hypothetical  organoli- 
thium  molecule  C2H2Li2  [9].  The  latter  hypothesis  has  been 
confirmed  by  semi-empirical  computer  simulations  on  model 
hydrocarbons;  in  particular,  we  find  no  evidence  for  stable 
epitaxial  Li2  molecules  [7]  nor  for  large-scale  occupancy  of 
first-neighbor  hexagons  by  intercalated  Li. 

Molecular  vibrational  spectroscopy  based  on  inelastic  neu¬ 
tron  scattering  offers  several  advantages  relative  to  the  more 
familiar  infrared  (IR)  and  Raman  spectroscopies,  in  partic¬ 
ular  very  high  sensitivity  to  modes  involving  H  motion  (due 
to  the  huge  incoherent  thermal  neutron  cross  section  of  the 
proton),  and  the  absence  of  symmetry  selection  rules  such 
that  all  modes  are  detectable.  We  used  the  reactor-based  filter- 


analyzer  neutron  spectrometer  (FANS)  at  NIST  [10]  to 
measure  the  vibrational  densities  of  states  of  graphite  and  the 
two  prototype  epoxy-based  carbons  described  above.  The  top 
panel  of  Fig.  4  shows  spectra  for  ENR650  (closed  circles) 
and  ENR1000  (closed  squares)  at  16  K.  Intense  features  in 
the  former  at  1 10  and  126  meV  are  assigned  to  out-of-plane 
vibrations  of  C-H  species,  with  neighboring  H’s  vibrating 
in-phase  and  out-of-phase,  respectively.  The  broad  maximum 
at  150  meV  originates  from  C-H  in-plane  vibrations.  These 
are  also  observable  in  ENR 1000  but  with  much  weaker  inten¬ 
sities,  and  are  ( sensibly )  absent  in  graphite.  The  lower  panel 
of  Fig.  4  shows  calculated  spectra  based  on  semi-empirical 
potentials  for  a  planar  7-hexagon  graphene  fragment  with 
each  edge  C  terminated  by  a  single  proton  ( coronene,  solid 
curve )  or  2  protons  ( dashed  curve ) .  The  experimental  results 
are  clearly  in  better  agreement  with  the  former,  confirming 
that  we  have  single  H  saturation  of  dangling  bonds  at  edge 
carbons.  Disordered  carbons  obtained  by  low-temperature 
pyrolysis  of  organic  precursors  can  therefore  be  reasonably 
modelled  as  a  distribution  of  large  polyaromatic  hydro¬ 
carbons. 

Addition  of  Li  to  these  C-H  entities  would  induce  a  down¬ 
shift  of  the  mode  frequencies  relative  to  the  C-H2  simulation 
since  Li  is  heavier  than  H.  We  modelled  the  stability  of  C- 
H-Li  species  using  semi-empirical  quantum  chemical  meth¬ 
ods,  validated  for  small  systems  by  ab  initio  LDA  and 
quantum  MD  [  1  ] .  An  example  is  shown  in  Fig.  5.  The  hex¬ 
agons  in  pyrene  C16HH)  are  terminated  by  either  2  or  3  hydro¬ 
gens  ( ‘2H’  and  ‘3H’).  Regardless  of  starting  configuration, 
a  single  Li  always  relaxes  to  a  position  approximately  cen- 
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Fig  4.  Top:  (•)  vibrational  densities  of  states  for  ENR650  and  (■) 
ENR  1000.  Bottom,  calculated  spectra  for  a  7-hexagon  graphene  fragment 
with  edge  carbons  terminated  by  ( - )  1 H  or  {  )  2H 
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Fig.  5.  Relaxed  structures  obtained  by  adding  2  Li's  to  C)f,H|((  on  opposite 
sides  of  the  molecule  centered  over  (top)  3H  and  (bottom)  2H  rings. 


tered  above  2H  (88.7  kcal/mol)  or  3H  (85.3  kcal/mol). 
Two  Li' s  initially  above  and  below  a  3H  remain  stable  ( Fig.  5 
top),  while  when  started  above  and  below  a  2H,  they  move 
off  the  interstitial  sites  and  bond  to  the  peripheral  H-termi- 
nated  carbons  (Fig.  5  bottom).  This  requires  sp3  hybridiza- 

c 

tion  of  both  carbons  (C-C  bond  length  =1.52  A),  and 
destroys  the  conjugation  of  that  particular  aromatic  ring  ( 1 ,2- 
bis-localization) .  Energy  differences  between  interstitial  and 
(C,H) -bonded  Li’s  are  insignificant.  We  therefore  expect 
both  to  exist  in  real  materials  at  similar  potentials,  strongly 
suggesting  that  residual  hydrogen  in  low-temperature  pyro- 
lyzed  amorphous  carbons  provides  a  second  mechanism  for 
Li  uptake  which  operates  in  parallel  with  intercalation  of  Li 
into  interstitial  sites.  Similar  results  are  found  for  larger  model 
hydrocarbons.  By  adding  Li’s  sequentially  to  a  large  poly¬ 
aromatic  hydrocarbon  and  fully  relaxing  the  intermediate 
structures,  we  find  that  Li’s  first  occupy  all  available  H- 
terminated  edge  sites  which  meet  the  criteria  described  above, 
then  go  to  second-neighbor  intersititial  sites  as  in  bulk  LiC6, 
and  finally  begin  to  form  isolated  Li  clusters.  At  no  point  do 
we  find  stable  epitaxial  Li2  molecules  [7]  nor  Li  multilayers 
[8]. 

4.  Stabilizing  LiC2 

High-pressure  reactions  of  Li  with  graphite  yield  LiC2, 
three  times  as  dense  in  Li  as  the  usual  LiC6  [11].  LiC2 
decomposes  rapidly  but  incompletely  into  LiC6-(-Li  upon 
releasing  pressure.  We  observe  a  dramatic  enhancement  in 
the  ( meta) stability  when  the  graphite  contains  0.5  at.%  sub¬ 
stitutional  boron  |2].  Average  compositions  derived  from 
X-ray  (OOL)  intensities,  determined  several  months  after 


Time  after  release  of  pressure  (days) 


Fig.  6.  Time  dependence  of  x  in  high  pressure-synthesized  LiC\,  measured 
in  days  after  releasing  the  pressure:  (O)  sample  prepared  from  ordinary 
HOPG  and  (  >)  sample  prepared  from  HOPG  containing  0.5  at.%  boron 
Horizontal  arrow  indicates  LiC6  for  reference. 


releasing  the  pressure,  are  LiC2  2  and  LiC3  8  with  and  without 
boron  respectively,  as  shown  in  Fig.  6.  The  origin  of  this 
phenomemon  is  not  known.  Boron  is  an  acceptor  in  graphite, 
depressing  the  Fermi  energy  and  increasing  the  electron  affin¬ 
ity.  But  it  is  hard  to  believe  that  1  boron  per  200  carbons 
could  have  such  a  marked  effect  on  the  Li  uptake  by  an 
electronic  mechanism  [12].  We  speculate  as  follows.  The 
metastable  composition  LiC2  2  with  boron  is  close  to  that  of 
Li 

27^72’  one  of  the  intermediate  decomposition  products 
[11].  Further  decomposition  requires  long-range  Li  diffu¬ 
sion.  Since  the  volume  per  boron  impurity  is  within  a  factor 
of  ~  2  of  this  unit  cell  volume,  we  suggest  that  boron  arrests 
the  decomposition  cascade  by  pinning  this  metastable  super¬ 
lattice  with  ~  1  boron  per  cell.  It  may  be  possible  to  exploit 
this  metastability  in  the  design  of  metal-free  high  capacity 
primary  batteries. 
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